Pulsar timing and spacetime curvature 
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ABSTRACT 

We analyze the effect of weak field gravitational waves on the timing of pulsars, with 
particular attention to gauge invariance, that is, to the effects that are independent of 
the choice of coordinates. We find: (i) the Doppler shift cannot be separated into gauge 
invariant gravitational wave and kinetic contributions; (ii) a gauge invariant separation 
can be made for the time derivative of the Doppler shift in which the gravitational wave 
contribution is directly related to the Riemann tensor, and the kinetic contribution is 
that for special relativity; (iii) the gauge dependent effects in the Doppler shift play 
no role in the program of gravitational wave detection via pulsar timing. The direct 
connection shown between pulsar timing and the Riemann tensor of the gravitational 
waves will be of importance in discussions of gravitational waves from alternative (non- 
Einsteinian) theories of gravitation. 

Subject headings: pulsars:general, gravitational waves 



X 



1. Introduction 

The poss ibility of dete cting gravitational wa ves through pulsar timing, first suggested inde- 
pendently by ISazhinl (|1978l ) and iDetweilerl (|1979l ) , is of increasing interest as part of the birth of 



gravitational wa ve astronomy, and the details of this technique continue to be advanced by sev- 
eral researchers (Hellings fc Downs 19831; backer h Hellings 1986 ; Hellings 1986 ; Kaspi et al. 1994 ; 

2008). The basis of this detection technique is 



Lommen fc Backer 



2001 



Jenet et al 



the eff ect that a gray i tational wave has on the arrival times of pulsar signals. In the analysis of this 



effect (IHellings 



1981 



Backer &; Hellings! Il986l ). the gravitational wave has been described in terms 



of perturbations h^ v = g^ v — of the spacetime metric g^ u . This description of gravitational 
waves is analogous to describing electromagnetic waves using the potentials {<!>, A}. In particular, 
coordinate transformations are analogous to the gauge transformations of electromagnetic theory, 
and - at least in linearized general relativity - are also called gauge transformations. 
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If electromagnetic fields are to be described in terms of potentials, then one must choose a 
particular gauge for the description. Similarly in relativistic gravity, coordinate conditions ("gauge 
fixing") must be added. The need for a gauge choice gives rise to the possibility of confusing a gauge 
effect with a physical effect. This confusion, for relativistic gravity, kept it a controversi al question 



for m any years whether gravitational waves were physical, or were coordinate waves (jKennefick 



2007). A further source of confusion has been the distinction between the Doppler shift due to 
gravitational waves, and the Doppler shift due to relative motion of the emitting pulsar and the 
receiving telescope. This has the potential to be particularly confusing, since the relative motion 
is a manifestation of spacetime curvature, as is the gravitational wave. 

Our purpose here is to clarify the description of pulsar timing, by analyzing pulsar timing 
with particular attention to what aspects of the pulsar Doppler shift are gauge dependent and 
what aspects are physical. What we will show is that pulsar timing can be analyzed in a way 
in which there is an incontrovertible contribution due to gravitational waves, and a contribution 
due to relative motion of emitter and receiver. The result is couched completely in terms of 
quantities that are intuitively appealing, as well as mathematically well defined. The gravitational 
wave contribution will be based on the Riemann tensor which is gauge invariant very much as 
are the electric and magnetic fields in classical electromagnetism. We will see, moreover, that in 
the literature there are claims about Doppler shift effects that are not gauge invariant, but that 
these statements refer to intermediate steps in the analysis of pulsar timing and do not affect the 
overall program of pulsar timing. The technique of gravitational wave detection by pulsar timing is 
physically valid; it is cleanly distinguishable both from coordinate effects and from source/receiver 
motions. 

Besides clarifying what is and is not physically meaningful in pulsar timing, this paper helps to 
establish a firm background for studying the nonstandard polarization modes of alternative theories 
of gravitation, modes for which pulsar-timing detection may be particularly useful. Such modes 
are most easily described in terms of components of the Riemann tensor, not in terms of metric 
perturbations. 

The paper is organized with Sec. [2] as its mathematical heart. In that section we derive the 
expressions that form the basis of our analysis, and that show clearly what is and what is not gauge 
invariant in pulsar timing. To help focus on the distinction between kinematic and gravitational 
effects, in Sec. [2] we assume that the motions of the emitter and the receiver of pulses are driven 
by nongravitational forces. We remedy this unphysical assumption in Sec. El where we add gravity 
as the source of the astrophysical motions of the emitter and receiver. In Sec. H] we discuss the 
implications of these results for the way in which pulsar-timing gravitational wave detection is to be 
carried out, and is to be viewed. In the Appendix we fill in some details of Sec. [2] that we thought 
would divert too much attention from the main point of that section. 

We have chosen to present these results with a minimum of unnecessary generality that would 
make the mathematics more elegant, but more obscure. Except as noted, we use the notation and 
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conventions of the text bv lMisner et al.l (|1973l ). 



2. Doppler shift and gauge invariance 

The Doppler shift from emitter to observer can be considered to have two sources: (i) general 
relativisitic: the effects of gravitational waves, and (ii) special relativistic: the relative motion of 
emitter and receiver. The second effect will be much larger than the first in astrophysical situations. 
The relative velociy of astrophysical bodies is on the order of v ~ 10~ 3 (in c = 1 units), while the 
characteristic magnitude of gravitational waves is smaller than 10 -20 , usually much smaller. In our 
development we keep special relatistic effects to all orders of v, but we will ignore effects of order 
v\h^ u \. The justification is that gravitational effects will be at the limit of detectibility; effects that 
are smaller by 10 -3 are not of immediate interest. Here and below we will use v to denote an 
astrophysical velocity, and we will repeatly use the fact that v\h^ u \ can be ignored. 

As pointed out in Sec. [H here we will make the artificial assumption that the emitter and 
receiver are point particles that are being driven in accelerated motions by nongravitational forces, 
such as rocket engines. In the current section, then, the only curvature of spacetime is due to 
gravitational waves. 

We will denote our emitter worldline as E, and receiver worldline as R. We choose Minkowski 
background coordinates t, x, y, z for the background so that they are appropriate for the approxi- 
mations just discussed. That is, in these coordinates the E and R worldlines are at rest aside from 
velocities of order v, and the metric perturbations /i™ due to gravitational waves are extremely 
small. The 4-momentum of a photon from emission to reception is written as = P°(l,n), so 
that n plays the role of a unit vector pointing (in the Minkowski background) from the emitter to 
the receiver. 

We now consider the following expression 

Dopp = - J [I htt,t + n 3 h tj<t + ±n J n fc h jkit J dX + - n k U k \ - |t7* - n k U k \ . (1) 

Here is the 4- velocity of of the emitter (E) and receiver (R), at the events of emission and 
reception; dX indicates integration along the photon worldline, with dt = dX and dx 3 = n^dX; Latin 
indices are spatial (referring to the x,y,z components of the coordinate basis). 

In the Appendix we show that the expression in Eq. ([1]) represents the Doppler shift of the 
photon, that is, the fractional difference by which the photon energy observed by the receiver is 
greater than that observed by the emitter. Here we focus on the gauge property of the expression, 
the changes induced in the expression by a coordinate transformation x^ ncw = x^ + £^ in which 
the gauge vector £ is of the order of the metric perturbations . 

The standard gauge transformations of the perturbations and of the components of the 4- 
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velocities are 

h^7 = v + w + u» now = eu» a - u a e, a ■ (2) 

We now note that the 4-velocities have the components = 8^ aside from corrections of order v. 
With terms of order v\£\ ignored we are left with 

U* 1 ncw = U!*-£%. (3) 

We also note that in Eq. ([1]) a gauge change in n would give terms of order v\£\ and s ° n ° 

gauge transformation of n is carried out. 

A straightforward calculation shows the gauge invariance of Dopp: 
5Dopp = Dopp ncw -Dopp = ~ J (e*,*,t + n'ttjt + n?t jM + n>n% k J) d\+ [& jt + n j ^ t ] R - + n'&t] R 

= -£(§- t + nP -^) (&.« + dA + + ^'^] R - [e« + ^,t] R = o • (4) 

(It is worth noting here that time derivatives of n can be ignored since these time derivatives are 
of order v and would be multiplied by terms of order |£|.) 

Since the Doppler shift, as defined, refers to an objective physical measurement, the fact that 
it is gauge invariant is simply a consistency requirement. The details of the gauge transformation, 
however, underscore an important point: neither the integral nor the 4-velocity contributions to 
Eq. (pQ) is separately gauge invariant. Thus, the temptation to identify the integral in Eq. <pQ) as 
the gravitational wave contribution and the 4-velocity terms as the kinetic contribution must be 
avoided, since that identification has no invariant meaning. 

To arrive at a more physically useful expression we take the derivate of Eq. (pQ) with respect to 
the coordinate time t. (In doing this we note again that dn/dt is of order v, and hence that time 
differentiation of n in the integral can be ignored.) The result of time differentiation is: 



dDopp 

dt 

+ 



J (\ h t t,tt + n j htj,u + 3 nPn k hjk,ttj d\ (5) 



dU l dW dn iTTi 

— 71,— -^-W 

dt 3 dt dt 



R 



dU l dW dn jT : 

— Tlj— -^-U 3 

dt J dt dt 



E 



A few comments here on the time differentiation are appropriate. The total time derivative of the 
integral should, in principle, include the change in the integral due to the change of the time of 
the end points of the integral. But this change involves v, and hence terms of order vlh^l which 
we ignore. The time changing endpoints, on the other hand, cannot be ignored in the 4-velocity 
terms, since these terms are not multipled by metric perturbations. The time derivatives d/dt 
in the 4-velocity terms are therefore understood to be the derivatives along the worldlines, i.e., 
d/dt = d t + v k d x k. 
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Next we consider the components of the 4-acceleration a M of the E and R worldlines 



7 



dU 3 
~~dT 

dU 1 



dU 3 , 
7 ~~jt~ + hj t ,t 



. =7 _ + (c , rr < 



7' 



dt 
dlP 
~dT 



2 h tt,j 



(6) 
(7) 



where we have ignored terms of order v\h^ u \ in the T term, and where 7, as usual, represents 
,2 When these equations are used in Eq. ([5]) we get 



1/vT 



ciDopp 
eft 



R 



n j n k h tt j k + n j n k h tjjt k 



\ n j n k h jkjtt 



d\ 



(8) 
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1 fl 


a* — Bja- 7 
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+ 
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3 dt _ 




dt 



n 3 n k R tjtk d\ + 



nM J 



7 



UjO? 



7 



+ 



-17, 



(in- 7 
~eft~ 



dn 3 
~dt 



(9) 



where Rtjtk represents the components of the Riemann tensor. 

The expression in Eq. (|9|) is, of course, gauge invariant, but unlike the gauge invariant expres- 
sion in Eq. ([1]), the individual contributions are now gauge invariant. The integrand contains only 
a projection of the manifestly gauge invariant Riemann tensor. The 4-velocity and 4-acceleration 
terms are all of first or higher order in v, so their gauge changes would be of order v \hav\, and hence 
ignorable. Unlike the the expression in Eq. ([1]), for the Doppler shift, the expression in Eq. Q for 
the time derivative of the Doppler shift contains contributions that have physical gauge-invariant 
meaning; the integral gives the effect of gravitational waves and the remaining terms give the special 
relativistic contributions due to acceleration (the 4-acceleration terms) and to the relative geometry 
of the worldlines (the dn/dt terms). 



3. Gravitationally driven orbits 

In Sec. [2] we assumed that the motions of the emitter and receiver were driven by nongravi- 
tational forces ("rocket engines"), so that the role of spacetime curvature lay solely in the pertur- 
bations hav identified as gravitational waves. In reality, of course, astrophysical orbital motions 
are driven by gravity. In this section we explain how to incorporate other aspects of gravity, in 
particular, orbital forces, into the considerations of Sec. El 

We imagine now that the emitter and receiver are on astrophysical bodies that are moving 
under the influence of other astrophysical bodies. We have already assumed in Sec. [2] that orbital 
velocities are small compared to c. This means that the gravitational interactions among all bodies 
are weak-field interactions. More precisely, for our system of gravitationally interacting bodies, the 
Newtonian potential $ due to one body, at the position of the other body, must be small compared 
to c 2 = 1. We now make the additional assumption that the gravitational field is weak (|$| -C 1) 
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everywhere along the photon path. This allows us to treat all non-gravitational-wave fie lds relevan t 



to the problem as being adequately described by a post-Newtonian (pN) approximation (|Willlll993l ). 
Though we may consider a higher order pN approximation, we do not consider orders high enough 
(order (v 7 )) for the orbiting bodies themselves to become sources of gravitational radiation. 

In a pN approximation, a background Minkowski-like coordinate system is used, and metric 



pertu rbations in terms of these potentials are required to have "post-Newtonian" character (jWill 



1993). This character requires, among other things, that the potentials be functions only of the 
coordinate separation of source point and field point. The equations determining the metric per- 
turbations are computed from Einstein's equations truncated to the order of the pN approximation 
being used. Gravitationally driven motions are then taken to be the geodesies of the pN metric, 
and may be interpreted as having accelerations in the fiat background. 

Due to the nature of the pN approach, a pN gauge vector would have to have pN character. 
This very strong constraint turn s out to le ave almost no gauge freedom. Moreover, there is a 



standard pN coordinate system (jWilll Il993l ) in which even the small gauge freedom is removed. 



When this gauge is adopted, the question of pN gauge transformations ceases to exist. 

Gauge transformations can still be made, of course. The h^ v gravitational wave metric pertur- 
bations do not have pN character. So long as our gauge field £^ is not pN in character it is clearly 
separable from any choice of pN gauge, and does not affect the construction of pN coordinates. 
Again we note that this presumes that we ignore any contribution to the field of gravitational 
waves from the pN-modelled bodies themselves. 

With these considerations we can conclude that the gauge transformations of Sec. [2] can be 
repeated in the pN plus gravitational wave spacetime, with the accelerations taken as those (relative 
to the Minkowski background) coming from pN theory. One might worry that the gravitational wave 
gauge transformations in Sec. [2] were made relative to a flat background, while the pN spacetime 
is not flat. But this simply means that we are ignoring terms of order |3>||£|, which are analogous 
to the terms of order t> 2 |£| that we ignored in Sec. [2j 

Though we are primarily interested here in gauge questions, the matter of the pN metric 
perturbations raises an important separate issue. The photon's Doppler shift will be affected by 
the surface gravitational field of the emitter and receiver. In addition, if the the photon happens 
to pass close to a another astrophysical body the pN fields of that other body wi ll affect pulsa r 



time-of-arrivals by altering photon path lengths and through the Shapiro time delay (|Shapirdll966l ) . 
These effects can be separately evaluated (as part of a pN calculation) and added to the Doppler 
shift calculations of Sec. O Since the effects are small, they do not "mix" with the terms. 

We close this short section by pointing out that there are limits to the clean separation of 
kinematic and gravitational wave terms. If the gravitational fields are strong, if orbital velocities 
are relativistic, or gravitational potentials along the photon path are of order unity, then the 
gauge-invariance demonstration in Sec. [2] fails to hold. Indeed, it is intuitively appealing that in 
such a case it should not be possible to make a general distinction between small perturbations 
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of spacetime and motions of astrophysical bodies. It would sill be possible, of course, to make a 
practical distinction between the two if the gravitational waves had a significantly higher frequency 
than the time scale of the kinematics. 



4. Conclusions 



A physical measurement, such as the Doppler shift of a photon, is not changed by the changes 
in our mathematical choices, so the gauge invariance of the expression in Eq. (fTJ) is simply a check 
of consistency. What is subject to our mathematical choices is the interpretation of the terms in 
Eq. ([1]), in particular the interpretation of the integral term as the gravitational wave contribution. 
This interpretation is clearly incorrect since we could, for example, choose coordinates that makes 
the integral vanish for any photon, and put all of the Doppler shift into the kinetic term. 



Though integrals like that in Eq. ([!]) have appeared in the literature (|Detweilerlll979l ; iMashhoon 



1982 



Bertotti et al.lll983l ) in discussions of pulsar timing, in practice, this causes no real difficulty. 



The Doppler shift can be understood to be the time integral of the gauge invariant expression 
in Eq. ([5]). That time integral will contain an integration constant that cannot be meaningfully 
separated into gravity and kinematics. The observed phase of the arriving pulses, the raw data of 
pulsar timing, will be the next time integral, and will contain a term A + Bt, where A and B are 
such integration constants. These integration-constant terms play no role in the actual analysis of 
timing residuals. The program of gravity wave detection by pulsar timing, therefore, has a solid 
physical foundation based on the Riemann tensor. 

This relationship of timing residuals to the Riemann tensor has a useful secondary benefit. As 
pointed out in the introduction, pulsar timing has the potenti al to be a particu larly sensitive probe 
of non-Einsteinian polarization modes of gravitational waves (|Lee et al.l 12008) . The description of 
these modes uses the components of the Riemann tensor. The expression in Eq. ([5]) gives a clear 
and unequivocal description of how these nonstandard gravitational waves affect pulsar timing 
residuals. 
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Fig. 1. — Construction of the timelike geodesic worldlines E g and R g , and of the Minkowski-like 
coordinates. 

Appendix: Derivation of the gauge invariant expression for the Doppler shift 

We now turn to a proof, needed in Sec. [21 that the expression in Eq. ([1]) is the Doppler shift. 
Since we have shown that the expression in Eq. ([1]) is gauge invariant, it suffices to show that it is 
equal to Doppler shift in any one gauge. 

We construct a convenient gauge for this proof as follows. We let e and r be the emission and 
reception events for a pulsar photon. We choose a timelike geodesic worldline E g through e to be 
almost tangent (i.e. , tangent to order v), at e, to the emitter worldline. Next we choose a timelike 
geodesic worldline R g through r in such a way that it is parallel to the geodesic worldline through 
e in the following sense. Near the emission event we construct a spacelike surface orthogonal to the 
emitter 4- velocity, and we construct this surface to have extrinsic curvature with a vanishing trace. 
We generate (much in the manner of Gaussian normal coordinates) a congruence of timelike geodesic 
worldlines normal to this surface. We assume that the spacetime curvature is small enough, and/or 
the emitter-receiver distance is small enough that this congruence fills spacetime, with no crossings, 
in the neighborhood of the reception event. We choose R g to be the curve in that congruence that 
goes through r. 

The next step in the construction is to generate null geodesies (photon worldlines) from E g to 
R g and to define our Minkowski coordinates t,x,y,z by the following steps: (i) the surface spanned 
by the null geodesies is taken to be a surface of constant x and y; (ii) z is set to zero along E g , and 
along R g we set z equal to another constant, the length of the spatial geodesic between E g and R g 
on the spatial hypersurface. (hi) z is taken to be an affine parameter along the null geodesies; (iv) t 
is taken to be proper time along E g ; (v) the coordinate t is propagated through the t, z surface by 
requiring t — z to be constant along the null geodesies, we note that a tangent to the null geodesies 
is given by 

^*L=const + ^zlt=const > (10) 
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in the coordinate system we have defined. 

We now note that the covariant t component of the photon 4-momentum satisfies the geodesic 
equation 

dP t _ dx 13 _l JD adxl 3 h /ill 



dX " dX 



and 



where 



d d 
dX~ dt 



z=const 



d_ 

dz 



t=const 



AP t = f P a dx t3 h a ^ t =P* [ ah tt , t + Htzjt + \h ZZit ) dX (12) 

JE JE 

(13) 

From this expression we have that the difference between the photon energies at r and e is 

f R 

[Energy measured on ~R g — Energy measured on E g ]/P t = — {^htt.t + htz,t + \h zz ,t) dX . (14) 

Je 

We now generalize this to the case in which the photon is not confined to the plane with x and y 
constant. Since the expression on the right in Eq. (114j) is already first order in the metric pertur- 
bations, we need to consider only a more general photon direction in the Minkowski background. 
If we denote by n the unit vector pointing in the spatial direction in which the photon moves, then 
the generalization of Eq. (fTH is 

[Energy measured on R g — Energy measured on E fl ]/P* = — (i^tt,t + n ^htj,t + \h zz $) dX . 

Je 

(15) 

To get the full expression for the Doppler shift we must consider the fractional energy changes 
from the geodesic worldlines to the observer worldlines. An observer with 4- velocity observes 
a photon with 4-momentum P^ to have energy —PfJJ^. By construction, our geodesic worldlines 
have components t/^ = {1,0,0,0}, Thus the energy observed by the receiver is P t [U t — n^Uj\ R , 
where P is the energy observed at the reception event by the geodesic observer. With this and the 
similar expression for the emission event we get Eq. ([1]). 
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